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Abstract 

Starting from the old idea that Fermi statistics for quarks play a fundamental 
role to explain some features of hadron structure, we study the modification of 
the scaling behaviour of parton distributions due to quantum statistical effects. 
In particular, by following an interesting formal analogy which holds between the 
Altarelli-Parisi evolution equations, in leading-log approximation, and a set of 
Boltzmann equations, we generalize the evolution equations to take into account 
Pauli exclusion principle and gluon induced emission. 
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Deep inelastic experiments seem to be an inexhaustible source of information on the 
hadronic structure and continue to considerably improve our understanding of strong 
interaction dynamics. A measurement of proton and neutron F 2 {x) structure function 
performed by the NMC Collaboration at CERN |] suggests a rather large SU (2) flavour 
breaking in the sea quark . In particular they have obtained a determination, at very 
small x, for the difference 

M*) = f -lF? P (y) - Ft{y)} = \C dy [«(„) + u(y) - d(y) - d{y)\ , (1) 
Jx y 6 Jx 

finding X G (0.004) = 0.227 ± 0.007 . Thus, by extrapolating down to x — 0, they have 
estimated 

J G (0) = 0.240 ± 0.016 . (2) 

This result represents a relevant violation of the Gottfried sum rule P| , which would 
predict X G (0) = 1/3. Moreover, from (|^) we get 

d - u = I dx \d(x) - u{x)\ ~ 0.14 . (3) 
Jo 

However, the inequality d > u was already argued many years ago by Field and Feyn- 
man [[| on pure statistical basis. They suggested that in the proton the production 
from gluon decays of im-pairs with respect to cM-pairs would be suppressed by Pauli 
principle because of the presence of two valence u quarks but of only one valence d 
quark. Assuming this point of view, the experimental result (0) naturally leads to the 
conclusion that quantum statistical effects play a sensible role in parton dynamics and 
that, in particular, parton distribution functions are affected by them. In this picture 
one may also easily account for the known dominance at high x of u-quarks over d- 
quarks, whose characteristic signature is the fast decreasing of the ratio F^x) / F^x) 
in this regime. 

In recent papers [Q this idea has been developed, succeeding in making reasonable 
assumptions for various polarized parton distributions in terms of unpolarized ones, 
explaining the observed violation of Ellis- Jaffe sum rule 0, and giving a possible 
solution to the spin crisis problem [^]. Finally, using Fermi-Dirac or Bose-Einstein 
inspired form for parton distribution functions, a rather good agreement has been 
obtained with the experimental data on structure functions ||. 

The aim of this letter is to study the role of quantum statistical effects, namely 
Pauli exclusion principle and induced gluon emission, in the Q 2 -evolution of parton 
distributions. We will show that there is a quite close analogy between the well-known 
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Altarelli-Parisi (A-P) evolution equations M, in leading- log approximation, and a set 
of Boltzmann equations written for a dilute system of particles. This analogy will guide 
us in finding the generalized scaling when also quantum statistics have been taken into 
account. 

As well-known, the logarithmic dependence on Q 2 of the parton distribution mo- 
menta, predicted in the framework of perturbative QCD, has a simple and beautiful 
interpretation in terms of evolution equations for parton distribution functions. At 
leading-log level, the A-P equations can be written in the following way 

l-M^f^PBivApJ*-) , (4) 

dt 2tt Jx y V \yj 

where t = ln(Q 2 //z 2 ), fi is some renormalization scale and pA{x,t) denote the parton 
distribution functions (A, B =quarks, antiquarks and gluons). By defining 
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In 



a(/j, 2 ) 
a{t) 



(5) 



with b = (33 — 2nf)/ (12tc) (rif is the number of flavours), Eq. (|4]) becomes 



4~Pa(x,t) = f —Y,PB(y,r)P AB (-] . (6) 
dr Jx y z \yJ 

Note that the dependence on r of r.h.s. of (|) comes only through ps(y, t). 
In Eqs. (H) and (|[), PAB^x/y) stand for the splitting functions evaluated by using 
standard equivalent parton method. They correspond to the probability for the ele- 
mentary three-body processes to occur in which a parton with momentum fraction x 
is produced by a parton with higher fraction y = x/z. 

The simple microscopical interpretation of Eq. ([]) is, as well-known, that the 
r dependence of pa distributions is induced by these processes considering them as 
occurring in the vacuum. 

However, it is physically reasonable to imagine that this picture has to be modified for 
sufficiently low x; in this regime the nucleons are filled with a large number of quark- 
antiquark pairs and gluons (the sea) and thus, to take into account in the correct 
way the presence of this large number of partons, the decays A —>■ B + C should 
be considered in presence of a surrounding plasma of both Fermi and Bose particles. 
Corrections induced by quantum statistical effects to the scaling behaviour dictated by 
(0) are therefore generally present, and in particular we expect that: 
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a) Pauli blocking will suppress the production of quarks and antiquarks with fraction 
x corresponding to filled levels; 

b) the gluon emission probability through bremsstrahlung processes, considered in 
the standard picture leading to A-P equations, will be enhanced by the contri- 
bution of induced-emission in presence of a rather relevant number of gluons in 
the sea. 

These effects would favour the production of gluon-quark pairs with larger values of x 
for the quarks and a smaller one for the gluon. Moreover the gluon conversion processes 
in q — q pairs are expected to be reduced. 

In statistical mechanics all these effects are simply included by multiplying the 
amplitudes modulus squared of the relevant processes by the factors 1 — / or 1 + / for 
each Fermi or Bose particle in the final state, with / denoting the particle distribution 
functions without any level-density factor. In equilibrium conditions these / reach the 
standard stationary Fermi-Dirac or Bose-Einstein form, while in general they depend 
on time. Thus, it is reasonable to expect that similar factors should be introduced in 
the A-P equations. 

To go further on this point we will show, as already mentioned, that A-P evolu- 
tion equations can be formally viewed as Boltzmann transport equations for parton 
distributions when they satisfy classical statistics (dilute system). 

As well-known, the Boltzmann set of equations describes the evolution to equilib- 
rium states of systems composed by many particles of several species (i specie-index 
i = 1, .., n) mutually interacting fllPl . Assuming, for the following applications, monodi- 



mensional dynamics for particles we can define the numerical distribution functions as 

m(e,t)= gi(e)fi(e,t) , (7) 

with e denoting the energy, fi(e,t) the statistical functions (they recover the usual 
Bose/Einstein or Fermi/Dirac at the thermal equilibrium), and <?i(e) the level-densities 
(weights) corresponding to e. These last quantities should be fixed from the beginning, 
by studying the hamiltonian of the total system. From (0) follows the expression for 
the total number of i-particles 



Ni(t) = Jde 9l (e) fl (e,t) . (8) 
By using Eq. ([?]), the Boltzmann equations can be cast in the following form 

£^ = a[f,g]=C+[f,g]-Cr[f,g] • = l,....,n , (9) 
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where f = (/i, /„), g = (gi, ■■■,g n ), C is the Liouville operator, and Cj[f, g] is the so 
called collisional integral for the i-th particle specie. The latter is given by a thermal 
average of all possible processes which change the density of the i-th specie. Notice that 
in Eq. @ we have defined Cf[i, g] and C^[f, g] as the contributions corresponding to 
the interaction processes which create or destroy the i-th particle specie respectively. 
For simple three body processes A — ► B + C, B — > A + C, if we are interested in 
describing the modification of B population, the corresponding terms in (7b [f, g] are 
the following 

C+[f,g]-C B [f,g] = 2vr J Jde A de c {\M(A^B + C)\ 2 5(e A -e B -e c ) 
x n A (e A ,t) g B (e B ) [l±f B (e B ,t)] g c {^c) [1 ± fci^c, *)]} 
-2tt J J de A de c {\M(B^A + C)\ 2 6(e B -e A -e c ) 
x n B (e B ,t) g A (e A ) [l±f A (e A ,t)] g c {e c ) [1 ± f c (e c , t)]} 

(10) 

where f = (f A (e A ,t),f B (e Bl t)J c (ec,t)), g = (g A (e A ), g B (e B ), g c (e c )), \M\ 2 are the 
squared moduli of transition amplitudes and the sign in the final state factors is pos- 
itive/negative depending on the bosonic/fermionic nature of particles. In the limit of 
very small fa one has (1 ± fi) ~ 1, and assuming free particle states (g^ = const) the 
collisional term for very dilute systems is recovered. 

Assuming rif different flavours for quarks (j = l,...,nf) with elicity states (A = 
+,— ), we can rewrite the set of equations @) for polarized quarks (qj\), antiquarks 
((jjx), and not-polarized gluons (G) distribution functions in terms of two-dimensional 
integrals 

d f 1 f 1 1 

■^Qjx{x,r) = / / dydz 8(x-yz) [P qq (z)q jX (y,T) + -P qG (z)G(y,T)} , (11) 

GST JO JO 2 

d r 1 f l 1 

—q jX (x,r) = dydz 5(x -yz) [P qq (z)q jX (y,r) + -P qG (z)G(y,T)} ,(12) 

d r 1 r 1 ( 

—G(x,r) = y J dydz8(x — yz) ^P GG (z)G(y,r) 

n f 

+ E E PcM 

j=l A=+- 

We have assumed for polarized gluon distributions that G + (x, r) ~ G_(x, r) ~ G(x, r)/2: 
we will comment on this point later. Note that in the previous expressions the 



Qjx(y,T) + q jX (y,T) } 



(13) 
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integrating-variables y and z vary from to 1, i.e. to the maximum available en- 
ergy properly normalized. Starting from the first of above equations, Eq. (p]), we 
notice that if we formally regard r as a time parameter, r.h.s. of ( [TT] ) represents 
the C^ x collisional term of a Boltzmann equation, written for particles obeying to a 
mono dimensional dynamics. This equation is given in terms of numerical distribu- 
tion functions qj\(x,r), G(x,r) and of the probabilities for the elementary processes 
Pqq(z) and PgG(z) (the final particles are assumed to be free which means no presence 
of extra-g terms corresponding to them). It is worth-while pointing out that in the 
infinite-momentum frame, where the parton picture is well-defined, all transverse dy- 
namics can be safely neglected (it has been already integrated out, indeed) and thus 
the description is mono dimensional. The (^-function in (|ll]) is just the longitudinal 
momentum conservation in the three body interaction process. 

However, to complete the analogy between (pTTj) and the corresponding Boltzmann 
equation we still have the difficulty that no Cj x terms are included in (|TJ): these cor- 
respond to decays of the x-momentum quark in parton pairs and would only depend 
on qj\(x, t) if all statistical factors for final states are neglected. The splitting func- 
tions involved in this case are two: P qq (z) and Pcq(z), being respectively, the emission 
probability of a quark or a gluon with fraction y = zx. Notice that in C~ x terms the 
momentum conservation leads to a definition of z which is simply the inverse of the 
one in C^ x . Thus, the C~ x collisional integral would be 

Cj x [q(x, r), G(x, r)] = q jX (x, r) £ £ dydz 5 (x - |) [P qq (z) + P Gq (z)} , (14) 

with q = (gi+, .., q nf +, <li-, •-, <?«/-), and by integrating over y we finally get 

Cr x [q(x,T),G(x,r)] = q jX (x,r) ^ dz z \P qq {z) + P Gq {z)} = , (15) 



where in fll5j ) we have used the integral constraint following from momentum conser- 
vation in quark splitting. Therefore from (^) since the C~ x contribution vanishes, we 
deduce that A-P evolution equation for qj X distribution can be consistently viewed as 
a Boltzmann equation in which all statistical factors corresponding to final particle are 
neglected (dilute system). Similar considerations can be repeated for the antiquarks 



evolution equations (|12|). 

A completely analogous result also holds for gluon distribution, whose evolution 
is dictated by (p!3|). In this case, by writing the corresponding Boltzmann inspired 
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equation, one gets 



A 11 f 

—G{x,t)= [ dydzS{x-yz){P GG (z)G(y,r) + J2 £ P g^) 

Oil JO JO j ^ ^ _|_ 

-1 r l 



x(y, r) + q jX (y, t)] } - G(x, r) J J dy dz 6 (x - |) [P GG (z) + 2n f P qG (z)\ 
= £-(p G c(-Wr) + £ £ P Gq (-)[ qj x(y,r)+q jX (y,r)]Y (16) 



which is just the A-P equation for gluon distributions. The C G term still vanishes due 
to the momentum conservation constraint 

f dz z[P GG (z) + 2n f P qG (z)}=0 . (17) 

J 

Let us briefly summarize our results till this point: we have shown that A-P evolution 
equations can be formally regarded as a set of Boltzmann equations for parton distri- 
bution functions, in which the Liouville operator takes the simple form of derivative 
with respect to r scaling variable. It is worth to notice that the absence of any external 
force in the regime of high transferred Q 2 is of course compatible with this expression 
for C. The analogy holds under the hypothesis that quarks and gluons form a very 
dilute system in the nucleons, so the statistical factors for final particles in the inter- 
action processes can be neglected. Starting from this equivalence and urged from the 
idea that, instead, quantum statistics would play a role in parton dynamics, it is now 
easy to generalize A-P equations to a set of generalized scaling equations where Pauli 
exclusion principle and gluon stimulated emission processes can be taken into account 
in a simple way. 

To this aim, as in equation fllOD , one should introduce in the collisional integrals 
the (1 ± fi) factors, and thus, as long as the statistical effects are taken into account, 
the factorization of qj\, cjj\ and G as reported in Eq. (|7|) becomes necessary. 
In the same spirit of (0), we will write the quark, antiquark and gluon distributions as 

qj\(x,r) = g jX (x) f}(x,r) , (18) 

Qj\( x , r ) = 9j\( x ) J*(x,t) , (19) 
G(x,t) = g G (x)f G (x,T) , (20) 

where gj\(x), g~j\(x) and g G (x) are weight functions, whereas fHx,r), fHx,r) and 
f G (x, t) are purely statistical distributions. The explicit form for (^-functions, which 
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contains the infrared divergency at x = 0, should be fitted from experimental data, as 
in P|, or deduced from theoretical expected behaviour, like, for example, Regge theory. 
We stress that the factorized form (]TE|)-(|211|), in particular the hypothesis that the 
singular functions gj\, gj\ and go do not depend on r is compatible with predictions of 
both Regge theory and QCD for the behaviour of parton distributions at the end-point 
x = 0; as well-known in this regime one has 



Pa(x,Q 2 )~UQ 2 



)x 



-OLA 



(21) 



with a.A which does not depend on Q 2 , at least for large Q 2 [ |TT[ . 
Remarkably, a parameterization similar to fll8"D - (pO|) has been already successfully pro- 
posed on phenomenological basis in M, to fit all the available measurements on parton 
distributions at fixed Q 2 , assuming for them a thermal-equilibrium form 



A x~ a 



A x~ a 



exp i ) + 1 



cxp 



x 



+ 1 
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G(x) = — A x~ a 



exp 



x 

x-x G 

X 



-1 



(22) 
(23) 
(24) 



where Xj\, Xj\, and x& represent the thermodynamical potentials, and x plays the 
role of the temperature. It is worth-while to point out that, in the framework of 
a formal connection between the A-P evolution equations and Boltzmann equations, 
these results || are a straight consequence of the above analogy, which predicts thermal- 
equilibrium-like solutions for parton distribution functions for sufficiently high values 
ofQ 2 . 

Within the factorized expression ([18|)- (|2"0D the final state factors are written in the 
form 1 — f*, 1 — and 1 + fc for quarks, antiquarks and gluons respectively 

We are now able to introduce a set of generalized scaling equations for quarks and 
gluons. Here we will consider for simplicity the case in which the gluons are supposed 
to not have a significant net polarization in the nucleons with respect to the one carried 
by quarks. We will assume, therefore G+(x,t) ~ G-(x, r) ~ G(x,r)/2. It should be 
pointed out that this approximation is consistent with the results obtained in || and 
H, where it is argued that still Pauli principle plays an essential role to generate the 
polarization of the quark sea, in the same approach therefore adopted here. 
It is, instead, less satisfactory in the framework of the different interpretation of the 
violation of Ellis- Jaffe sum rule based on the axial- vector current anomaly ||12j|. This 
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latter case, in fact, would require a very large gluon polarization, i.e. AG = G + — G_ ~ 
3 -j-4. Notice however that, as shown in ||, gluons are expected to be more numerous 
than quarks, due to their Bose nature, so in any case one has AG/G « Aq/q, which 
supports our approximation. 

By helicity conservation at the quark-gluon vertex, it is easily seen that evolution 
equations for polarized quark distribution functions get the following form 



d^ 



Qj\{x,r) 



1 dz 

x Z 



X 



l + -f G (X --1 ,T 



+ -P qG (z)G[-,r 



X 



/r A (^(--i)>r 



- Qjx(x, r) z dz\ P qq (z) 1 - /. (xz, r) 



1 + -/ g (x(1-z),t) 



+ Pg 9 (*) 



1 + -fG(xz,r) 



f}(x(l-z),T) 



(25) 



The equations for antiquarks are easily obtained by the previous one by substituting 
Qjx QjX an d fj <-» /j A . Similarly for the gluon distribution G(x, r) one has 



±G(x,t) 
dr 



ld ^h GG (z)G(^r 



l + -f G (x,r) 



T 



+ E E 



l + -f G (x,r) 



X 



1 + ^faixz, r) 



l + -/ (x(l-z),r) 



+ ^E E {[l-//(^,r)][l-^- A (x(l-«),r) 



+ 



[l-j^r)] [l-/- A (x(l-z),r) 



}} 



(26) 



Several comments can be made on the above expressions. First of all we notice that 
the inverse decay processes, the ones contained in the C~ collisional integral, contribute 
to the scaling behaviour of parton distribution functions with terms quadratic in the 
Pa{x,t) at least. This is a consequence of our interpretation of A-P equations as 
transport equations. As already shown these terms vanish in the limit of a very dilute 
system. 
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The generalized equations predict also a different, more complicated, evolution for 
momenta. By taking Mellin transform of both sides of ( p5|) and (|26|) , in fact, one sees 
that the standard scaling behaviour should be corrected by terms quadratic and cubic 
in distribution functions, which are not simply products of momenta of quarks and 
gluon densities. 

Finally, as for the standard A-P equations, the scaling behaviour for unpolarized 
quark distributions can be obtained by simply considering the sum qj(x, r) = qj + (x, r) + 
qj-(x, t) (the same holds for antiquarks). Notice, however, that since the introduction 
of final state statistical factors spoils the linearity of the equations, the evolution of 
qj (x, t) will depend on both the polarized distribution functions and not simply on 
their sum. 

We are now at some concluding remarks. We have stressed the point that, as 
some experimental results suggest, the Fermi or Bose nature of partons could sensibly 
manifest itself in observable quantities in deep inelastic scattering on nucleons. This 
idea already successfully applied in || and [§] mostly motivates our paper. In particular 
it seems to us quite natural that quantum statistics may modify the scaling behaviour 
of parton distribution functions for rather small x and high Q 2 ; in this region the sea 
becomes dominant and thus bremsstrahlung processes, responsible at leading-log level 
for scaling breaking, are likely supposed to occur in presence of a gas of partons. In 
this case Pauli blocking and gluon stimulated emission play, in general, a sensible role 
in parton distribution evolution. 

We have introduced both this statistical effects obtaining the generalized scaling, 
starting from the observation that a quite close and intriguing analogy seems to hold 
between A-P equations and a set of Boltzmann transport equations for partons. In 
our approach to nucleons as statistical systems, this fact is of course welcome and 
expected. Pursuing this formal analogy it seems to us very fascinating the fact that 
the scale variable Q 2 can be interpreted in some sense as a time parameter. The 
physical significance of this point, if any, should be deeper understood. It is also rather 
interesting to stress that, from this point of view, one would naturally expect, in the 
spirit of Boltzmann H theorem, that the normalized parton distributions fj\(x,r), 
fjx(x,r) and fc(x,r) should approach stationary Fermi and Bose expressions as Q 2 
increases. Remarkably, these conclusions seem to agree with the phenomenological 
results obtained in ||. 

We are much indebted with Prof. Franco Buccella, which we are pleased to thank 
for indefatigably encouraging this work and for his valuable comments. We would like 
also to thank Prof. E. W. Kolb, Dr. V. A. Bedniakov and Dr. S. G. Kovalenko for 
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